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ABSTRACT 


A  method  was  developed  for  forecasting  swell  using  a 
spectral  wave-generation  model  based  on  a  multiple  point- 
source  concept  of  swell  origin.   The  multiple  point-source 
concept  considers  that  the  peak-energy  swell  emanating  from 
a  moving  cyclonic  storm  can  be  considered  to  have  been  pro- 
duced at  one  or  more  space- time  point  sources  in  the  storm 
by  the  impulsive  introduction  of  energy  into  the  sea.   The 
method  was  tested  on  five  North  Pacific  storms  generating 
swell  recorded  at  Monterey,  California.   Predicted  swell  QfycicM^10' 
hej^ghts,  which  were  made  for  two  storms,  were  significantly 
lower  than  the  observed  heights.   The  time_  of  occurrence  of 
the  predicted  peak  height  agreed  with  that  observed  for  the 
swell  from  one  storm,  but  differed  by  about  ten  hours  for 
the  other.   Predictions  of  the  dominant  swell  period  jwere 
accurate  to  within  about  one  second  over  the  entire  range 
of  observed  periods  for  all  storms. 
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I.   OBJECTIVE  OF  THE  STUDY 

The  objective  of  this  study  is  to  develop  a  method  for 
forecasting  swell  at  a  single  observation  site  based  on  a 
multiple  point-source  concept  of  swell  origin  and  using  a 
spectral  wave  generation  model. 

The  multiple  point-source  concept  of  swell  origin  was 
proposed  to  this  investigator  by  Professor  Warren  C. 
Thompson  of  the  Department  of  Oceanography,  Naval  Postgrad- 
uate School.   The  concept  involves  the  assumption  that  the 
peak-energy  swell  emanating  from  a  moving  cyclonic  storm  can 
be  considered  to  have  been  produced  at  one  or  more  space- 
time  point  sources  in  the  storm  by  the  impulsive  introduc- 
tion of  energy  into  the  sea. 

The  swell  generation  model  proposed  is  a  spectral  model 
utilizing  the  Pierson-Moskowitz  (1964)  spectrum  and  assumes 
that  the  spectrum  of  peak-energy  waves  present  at  a  point- 
source  can  be  computed  from  the  surface  wind  speed  at  that 
point  and  the  speed  of  the  storm. 

The  products  of  the  method  are  a  forecast  of  the  domi- 
nant swell  period  and  significant  swell  height  with  time  at 
the  forecast  station. 

The  method  was  tested  on  five  North  Pacific  storms  and 
the  forecast  products  were  compared  with  the  observed  swell 
recorded  at  Monterey,  California. 


II.   DEVELOPMENT  OF  THE  METHOD 

A.   BASIC  APPROACH 

In  a  notable  study  of  swell  propagated  over  long  distances 
and  recorded  offshore  from  San  Clemente  Island,  California, 
Munk,  et  al.  (1963)  observed  slanting  ridges  in  the  energy- 
density  topography  drawn  on  a  plot  of  frequency  versus  time. 
Each  of  these  ridge  lines  described  the  frequency- time  dis- 
tribution of  the  maximum  energy  in  a  swell  train  arriving  at 
the  recorder  site  from  a  given  storm.   It  was  shown  that  a 
ridge  line  can  be  accounted  for  in  terms  of  classical  wave 
theory  so~  long  as  the  energy  spectrum  is  assumed  to  have  been 
generated  at  a  point-source  in  space  and  time.   The  slope  and 
zero-frequency  intercept  of  a  particular  slanting  ridge  line 
was  then  used  to  compute  the  effective  origin  time  and  dis- 
tance from  the  station  to  the  source  of  the  peak  energy 
swell. 

The  swell  prediction  model  developed  in  the  present  study 
effectively  employs  the  reverse  of  the  above  procedure.   A 
"point-source"  is  identified  in  the  wind  field  of  a  moving 
cyclonic  storm  on  each  of  a  succession  of  six-hourly  weather 
maps  covering  the  life  of  the  storm.   Each  point-source 
represents  a  potential  point  of  origin  of  the  dominant  wave 
energy  emerging  from  the  storm  and  is  considered  to  be  lo- 
cated at  that  point  where  the  peak  winds  are  generating  waves 
directly  toward  the  distant  station  for  which  the  prediction 


is  to  be  made.   A  wave  spectrum  is  computed  at  each  point- 
source;  the  characteristics  of  each  spectrum  are  determined 
from  the  surface  wind  speed  at  that  point  and  the  speed  of 
the  storm  toward  the  station.   FarV)  pprinH  mmpnnent  pf  t.hft 
spectrum  is  propagated  toward  the  station  at  a  group  veloc- 
ity which  is  a  function  of  its  period.   The  dispersive  ar- 


rival at  the  station  of  the  waves  propagating  from  a  single 
point-source  is  represented  by  a  plot  of  period  versus  time 
of  arrival. 

The  energy  in  each  spectral  component  of  the  swell  train 
represented  by  a  particular  period-time  curve  is  calculated 
by  modifying  the  energy  spectrum  at  the  point-source  for  the 
effects  of  angular  spreading  and  hiph-frequencv  attenuation.,^ 
By  graphing  the  energy-density  values  associated  with  the 
various  period  components   as  a  function  of  time  of  arrival 
at  the  station,  an  energy  density-time  plot  is  obtained. 
This  plot  is  termed  a  "propagated  spectrum".   The  envelope 
of  the  family  of  propagated  spectra  obtained  from  the  series 
of  point- sources  during  the  life  of  a  storm  represents 
the  energy  density  of  the  spectral  peaks  at  each  arrival 
time  at  the  observation  station. 

The  total  energy  in  the  swell  at  any  given  arrival  time 
is  estimated  by  computing  the  area  under  the  plot  of  fre- 
quency versus  energy  density  values  obtained  from  a  time- 
cut  of  the  propagated  spectra.   The  predicted  swell  heights 
are  assumed  to  be  statistically  related  to  these  coarse 
energy  calculations.   The  time  distribution  of  dominant 


swell  periods  is  obtained  from  the  time  distribution  of  the 
frequency  components  containing  the  maximum  energy  density 
in  the  propagated  spectra. 

In  the  following  sections  these  procedures  will  be 
described  and  illustrated  using  synoptic  weather  events. 

■B.   GENERATION  OF  SEA  SPECTRA 

1 .   Selection  of  Sea-Level  Pressure  Charts 

The  weather  charts  used  in  this  study  were  six- 
hourly  sea-level  pressure  analyses  of  the  North  Pacific 
Ocean  produced  by  the  Fleet  Numerical  Weather  Central  (FNWC) , 
Monterey.   An  advantage  of  using  these  computer-produced 
analyses  is  the  objective  manner  in  which  the  isobaric  fields 
are  produced.   The  fields  analyzed  covered  the  entire  area  of 
the  North  Pacific  under  investigation. 

The  weather  charts  were  chosen  on  the  basis  of 
frequency- time  analyses  of  five  selected  swell  trains 
recorded  at  Monterey.   Each  of  these  swell  trains  had  been 
analyzed  for  height  and  period  of  the  dominant  waves.   The 
storm  producing  each  train  was  identified  on  the  weather 
maps  from  the  effective  origin  time  and  distance  derived 
from  the  observed  frequency- time  distribution  in  the  swell 
in  the  manner  described  by  Munk,  et  al.  (1963).   All  were 
North  Pacific  storms.   The  sea-level  pressure  charts  cover- 
ing the  important  part  of  the  history  of  each  storm  are 
listed  in  Table  1,  along  with  the  corresponding  swell  trains 
recorded  at  Monterey. 


Table  I:   STORM  AND  SWELL  DATA 

A.   North  Pacific  Storms  Studied 

Storm 
Designation    Surface  Pressure  Charts  Analyzed  (GMT) 

1  1200/15  Feb  to  1200/17  Feb  67 

2  1800/17  Feb  to  1800/20  Feb  67 

3  1200/12  Nov  to  1200/15  Nov  67 

4  0000/16  Nov  to  0600/19  Nov  67 

5  1200/12  Mar  to  0600/15  Mar  69 

B.   Associated  Swell  Trains  Recorded  at  Monterey 

Swell  Train 

Designation    Observed  Swell  Occurrence  at  Monterey  (GMT) 

1  1200/20  Feb  to  1600/23  Feb  67 

2  0800/23  Feb  to  0800  25  Feb  67 

3  1200/17  Nov  to  2200/20  Nov  67 

4  0400/21  Nov  to  0600/23  Nov  67 

5  2000/18  Mar  to  0800/20  Mar  69 
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2 .   Determination  of  Surface  Winds 

By  measuring  the  spacing  between  isobars  on  sea-level 
pressure  charts  the  average  geostrophic  wind  speed  was  cal- 
culated from:  K     ap 

9    ~        2  9.  P   sin  $      An 

where   V   =  geostrophic  wind  speed  in  knots 

Ap  =  pressure  differential  in  millibars 
An  =  isobar  spacing  in  degrees  latitude 

6    =  latitude  in  degrees 

-  3       3 
p   =  density  of  air   (1.26  x  10    gm/cm  ) 

K  =  units  conversion  factor  (1.725  x  10   ) 

ft  =  earth's  angular  speed  of  rotation  (7.29  x  10 
radians/sec) 

In  this  study,  because  of  the  relatively  small  scale 
of  the  weather  maps  used  (1:60,000,000),  geostrophic  winds 
were  calculated  using  a  12-millibar  interval  (3  isobar  spac- 
ings) .   Geostrophic  winds  were  converted  to  surface  winds 
using  a  surface  to  geostrophic  wind-speed  ratio  of  0.8. 
FNWC  currently  uses  an  average  ratio  of  0.78  based  on 
Carstensen's  (1967)  report  showing  the  variation  in  the  ratio 
with  latitude  (Figure  1).   In  the  present  investigation,  the 
five  storms  studied  were  located  between  latitudes  31  and 
53  degrees  North;  accordingly,  a  ratio  of  0.8  was  considered 
to  be  a  more  representative  average  for  these  latitudes. 
The  surface  wind  direction  was  assumed  to  be  at  a  cross- 
isobar  angle  of  15  degrees;  this  agrees  with  current  Fleet 
Numerical  practice. 
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Figure  1:   Surface  Wind-Geostrophic  Wind  Ratio  (Vs/Vg)  as  a 
Function  of  Latitude  (after  Carstensen,  1967) 
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Figure  2:   Point-Source  Location  for  Surface  Pressure 
Chart  1200Z/15  Nov  1967 
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3 .  Location  of  the  Point-Source 

From  consideration  of  the  isobar  curvature  and 
surface  wind-speed  profile  characteristic  of  a  cyclonic 
storm,  it  is  assumed  that  waves  of  maximum  energy  propagat- 
ing toward  Monterey  must  be  generated  at  or  close  to  the 
point  where  the  maximum  surface  winds  are  directed  toward 
Monterey . 

In  order  to  identify  that  location  in  the  cyclonic 
wind  field  on  a  selected  weather  map, an  overlay  of  great- 
circle  arcs  from  Monterey  was  prepared  (Figure  2) .   By 
rotating  the  overlay  fifteen  degrees  clockwise  about  the 
storm  center  the  locus  of  points  in  the  storm  where  the 
great  circle  arcs  are  tangent  to  the  isobars  can  be  deter- 
mined; this  locus  represents  the  storm  radius  along  which 
are  found  those  surface  winds  generating  waves  directly 
toward  Monterey.   The  point  of  maximum  wind  speed  along  this 
radius  is  identified  by  measuring  successive  3  isobar 
(12-millibar)  spacings  and  locating  the  midpoint  of  the 
strongest  pressure  gradient.   Figure  2  illustrates  this 
technique  for  a  selected  surface  pressure  chart. 

4 .  Specification  of  the  Sea  Spectrum  at  the  Point-Source 
a.   Basic  Considerations 

Moskowitz  (1963)  analyzed  wave  records  and 
corresponding  wind  measurements  taken  by  ocean  weather  ships 
of  the  United  Kingdom  in  the  North  Atlantic  Ocean  and 
obtained  a  family  of  spectra  representing  fully-arisen  sea 
conditions  for  certain  wind  speeds  between  20  and  40  knots. 
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Pierson  and  Moskowitz  (1964)  later  used  the  data  to  test  the 
similarity  theory  of  Kitaigorodskii  who  proposed  that  all 
fully  developed  seas  should  have  the  same  non-dimensional 
spectral  form.   They  showed  that  within  the  limit  of  accu- 
racy of  the  reported  wind  speeds,  the  data  of  Moskowitz 
confirmed  the  theory  that  a  single  non-dimensional  power 
spectrum  might  be  used  to  represent  the  dimensional  spectra 
for  fully  arisen  seas  generated  by  any  wind  speed.   The  non- 
dimensional  spectrum  reduced  from  these  data  was  determined 
to  have  its  non-dimensional  peak  value  of  2.75  x  10    at  a 
non-dimensional  frequency  of  0.140  using  the  following  non- 
dimensionalizing  conversions: 

r=    IVs/g  (1) 


sir)  =  E(<)g3/vs5  (2) 


where   f  =  dimensional  frequency  in  hertz 
f  =  dimensionless  frequency 
E(0=  dimensional  energy  density  at  frequency  f 
..  S(f)  =  dimensionless  energy  density  at  frequency  T 
V$  =  surface  wind  speed  in  m/sec 
9  =  acceleration  of  gravity  (9.8055  m/sec  ) 

Bretschneider  (1963)  reported  that  the  general 
analytical  formula  for  a  one-dimensional  gravity  wave  spec 
trum  can  be  written: 


S(l)  =  K2  rm  e-(B2*~n)  (3) 
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where   S(l)   =  spectral  energy  density  for  frequency  f  ,  and 
where  it  can  be  shown  that: 


B   =  ™f 
2    n  TP 


K2=i;s|.p)e"'" 


f  =  frequency  of  the  spectral  peak 


In  a  later  work,  Bretschneider  (1963)  proposed  a 
spectral  form  for  wind- generated  seas  based  on  Equation  3  in 
which  m  =  5  and  n  =  4.   Pierson  and  Moskowitz  concluded  that 
for  the  purpose  of  forecasting  the  properties  of  the  larger 
waves  in  a  wind-generated  sea,  this  spectral  form  is  a  good 
representation. 

b.   Generation  of  Energy  Spectra 

In  this  investigation,  Bretschneider ' s  analytical 
formula  (Equation  3)  was  used  to  generate  Bretschneider ' s 
(1963)  spectral  form  approximation  of  the  Pierson-Moskowitz 
non-dimensional  spectrum  over  the  dimensionless  frequency 
range  0.10  to  0.30,  using  an  IBM  1620  computer.   Figure  3  is 
a  plot  of  this  spectrum. 

The  dimensionless  spectrum  was  used  to  obtain  a 
dimensional  power  spectrum  at  the  point-source  on  each 
weather  chart  using  the  wind  speed  value  calculated  for  the 
point-source  and  the  conversion  factors  of  Pierson  and 
Moskowitz  (Equations  1  and  2) . 
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Figure  3:   Dimensionless  Power  Spectrum 
(after  Bretschneider ,  1963) 
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The  spectrum  calculated  using  this  procedure 
represents  a  fully  arisen  sea  condition.   Since  the  sea  at 
the  point-source  on  most  weather  maps  was  clearly  not  fully 
arisen,  the  fully  arisen  spectrum  was  truncated  at  an  appro- 
priate cutoff  frequency,  as  proposed  by  Neumann  (1953),  to 
approximate  the  spectrum  of  the  non-fully  developed  sea  con- 
sidered to  be  present.   The  basis  for  the  procedure  devised 
follows  arguments  presented  by  Braunstein  (1970) . 

Braunstein  argued  that  fully  arisen  seas  should 
be  generated  when  the  velocity  of  the  fetch  equals  the  group 
velocity  of  the  peak  swell  being  generated,  since  such  a 


coincidence  produces  an  effectively  long  duration.   In  the 


present  investigation  a  group  velocity-fetch  velocity 
coincidence  was  used  to  effect  a  duration  limit  on  the  fully 
arisen  energy  spectrum.   If  one  travels  with  the  dominant 
waves  at  their  group  velocity,  the  energy  put  into  the  sea 
at  a  given  wind  speed  is  a  function  of  the  duration  over 
which  the  waves  remain  under  the  influence  of  the  generating 
wind  field.   In  the  case  where  the  waves  travel  in  the  direc- 
tion of  storm  movement,  as  is  the  circumstance  in  this  study, 
the  longest  effective  wind  duration  is  achieved  for  that 
frequency  component  whose  group  velocity  equals  the  speed 
of  movement  of  the  generating  source.   Lower  frequency  com- 
ponents, if  generated,  move  faster  than  the  generating  source 
and  move  ahead  of  the  area  of  peak  winds  before  they  acquire 
full  energy.   Higher  frequency  components  require  a  rela- 
tively short  wind  duration  for  achieving  saturation  and 
remain  in  that  state. 


In  this  study,  Neumann's  (1953)  model  of  a 
sharp  spectral  cutoff  was  assumed,  and  the  spectrum  gener- 
ated by  the  peak  winds  at  the  point-source  on  a  given 
weather  map  was  truncated  at  that  frequency  having  a  group 
velocity  equal  to  the  speed  of  the  point-source  averaged 
over  the  six-hour  period  just  preceding. 

This  method  has  some  obvious  deficiencies.   For 
very  large  storm  speeds,  it  tends  to  produce  spectra  which 
are  more  developed  than  would  actually  be  the  case,  and  for 
very  low-speed  storms  (including  the  stationary  storm  case) 
the  spe_C-tra  4mo_duced  are  less  developed  than  would  be 
expected.   As  will  be  shown,  however,  for  forecasting  the 
characteristics  of  the  dominant  swell  produced  in  a  moving 
cyclonic  storm,  the  method  produces  results  which  are 
reasonably  consistent  with  observations. 

The  problem  of  specifying  the  duration  associ- 
ated with  the  generation  of  the  maximum  seas  in  a  moving 
cyclonic  storm  of  given  size,  intensity,  and  speed  is  not 
simple  and  merits  additional  research. 

C.   SWELL  PROPAGATION 

1 .   Dispersion  of  Period  Components 

The  spectral  components  of  the  predicted  seas  are 
then  permitted  to  propagate  to  Monterey.   The  arrival  time 
at  Monterey  from  a  given  point-source  is  computed  for  each 
component  using  the  dispersion  relation: 

fa  =  f0  +  4  7rd/gT 
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where    *a  =  arrival  time  of  period  component  T 

*0  =  origin  time  at  the  point-source  (a  weather-map 

time) 
d  =  great  circle  distance  from  the  point-source  to 

Monterey 
T  =  period  in  seconds 

Figures  4  through  8  are  plots  of  swell  period  versus 
time  of  arrival  at  Monterey  for  each  of  the  five  storms 
studied.   For  each  storm  there  is  a  set  of  curves,  each  of 
which  represents  the  distribution  with  time  of  the  compon- 
ents of  the  wave  spectrum  arriving  at  Monterey  from  a  given 
point- source .   Each  curve  is  labeled  with  a  date-time  group 
indicating  the  map  time  of  the  particular  sea-level  pressure 
chart  containing  the  point-source  from  which  the  curve  was 
derived.   The  length  of  each  curve  is  governed  by  the  char- 
acteristics of  the  sea  spectrum  calculated  at  the  originat- 
ing point-source.   The  high  frequency  end  (that  component 
arriving  latest  in  time)  is  determined  by  introducing  a 
fixed  minimum  energy  density  criterion  into  the  dispersion 

program;  frequencies  in  the  sea  spectrum  having  energy  den- 

2 
sity  values  less  than  1.0  m  -sec  are  not  plotted.   The  low 

frequency  end  of  each  curve  represents  the  frequency  of 

spectral  cutoff  discussed  in  the  previous  section. 

2 .   Propagation  of  Spectral  Energy 

a.   General  Considerations 

In  addition  to  the  effect  of  dispersion,  the 

energy  in  the  wave  spectrum  generated  at  each  point-source 
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Figure  4:   Period-Time  Curves  for  Storm  1 
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Figure  5:   Period-Time  Curves  for. Storm  2 
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Figure  6:   Period-Time  Curves  for  Storm  3 
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Figure  7:   Period-Time  Curves  for  Storm  4 
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Figure  8:   Period-Time  Curves  for  Storm  5 
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is  diminished  by  the  effects  of  angular  spreading  and  high 
frequency  attenuation  during  propagation.   The  angular 
spreading  function  is  the  frequency- independent  Neumann 
function  (Figure  9) .   An  attenuation  function  is  used  to 
account  for  the  non-linear  processes  of  wave-wave  interac- 
tion and  wave  breaking;  it  results  in  dissipation  of  energy 
mainly  in  the  high-frequency  end  of  the  spectrum  with  little 
noticeable  effect  on  the  low  frequencies, 
b.   Angular  Spreading 

Although  on  each  weather  map  the  maximum  energy 
waves  being  produced  are  considered  to  have  originated  at  a 
point-source,  the  generating  area  has  a  finite  width  by 
virtue  of  the  fact  that  the  maximum  surface  wind  speed  com- 
puted at  the  point-source  is  obtained  by  averaging  the  wind 
over  a  12-millibar  isobar  interval.  The  effective  width  of 
the  generating  area  is  therefore  taken  to  be  equal  to  this 
12-millibar  spacing. 

An  angular  spreading  factor  is  calculated  using 
the  fetch  width,  the  great  circle  distance  from  the  point- 
source  to  Monterey,  and  the  trigonometric  relationships 
outlined  in  Pierson,  Neumann,  and  James  (1955),  as  illus- 
trated in  Figure  9.   Because  the  surface  wind  at  the  point- 
source  is  directed  toward  Monterey,  the  effective  fetch 
width  as  viewed  from  Monterey  has  a  fixed  angular  relation- 
ship to  the  great  circle  trajectory  from  the  point-source 
to  Monterey.   This  circumstance  simplifies  the  geometry  of 
angular  spreading  and,  along  with  the  adoption  of  a  fetch 
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width  defined  by  the  pressure  gradient,  makes  it  possible 
to  calculate  angular  spreading  factors  as  a  function  of 
distance  and  wind  velocity. 

For  practical  use,  a  table  of  angular  spreading 
factors  was  computed  for  wind  speeds  of  10  to  70  knots,  at 
increments  of  great  circle  distance  of  60  nautical  miles 
(for  distances  of  600  to  5000  NM) ,  and  at  latitude  incre- 
ments of  one  degree  (from  30  to  70  degrees  North) . 

In  the  present  study  the  swell  are  propagated 
spherically  along  great  circle  paths;  accordingly,  the  PNJ 
method,  which  assumes  straight- line  propagation,  was  modi- 
fied using  spherical  trignometric  relations.   A  comparison 
of  the  calculations  using  both  methods  indicated  that  the 
slight  increase  in  accuracy  obtained  using  spherical 
trigonometry  was  not  sufficient  for  the  propagation  ranges 
considered  (2500  to  5000  nautical  miles)  to  warrant  the  more 
complicated  spherical  computations. 

c.   High-Frequency  Attenuation 

Although  the  subject  of  wave-wave  interaction 
and  other  non- linear  processes  has  received  considerable 
treatment  in  the  recent  literature  (e.g.,  Hasselman,  (1963), 
Phillips,  (1963),  and  others),  no  analytical  representation 
has  been  derived  that  is  suitable  for  application  to  long- 
range  swell  prediction.   However,  Snodgrass,  et  al.  (1966) 
presented  empirical  data  on  the  attenuation  of  selected  fre- 
quencies observed  in  their  study  of  long-range  propagation 
of  ocean  swell  in  the  Pacific.   In  general,  they  found 
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attenuation  to  be  large  within  the  limits  of  the  wind  area 
of  the  generating  storm  and  negligible  beyond  the  storm. 
In  swell  that  had  travelled  a  considerable  distance  from  the 
storm,  they  also  found  that  below  a  frequency  of  0.06  hertz, 
attenuation  was  too  small  to  be  measured  and  above  0.0.8 
hertz,  the  wave  energy  was  masked  by  the  background  level. 
The  empirical  attenuation  data  given  by  Snodgrass,  et  al. 
were  logarithmic  attenuation  coefficients  reported  in  units 
of  decibels  per  latitude  degree  of  propagation  distance. 
For  the  range  of  frequencies  0.06  to  0.08  Hz,  these  data  fit 
an  energy  attenuation  function  of  the  form: 

-2ax 

e 

where   a.  =  modulus  of  amplitude  decay  in  deg-1  =  0.1151  /S 

/3  =  logarithmic  attenuation  coefficient  in  dB/degree 
x  =  propagation  distance  in  degrees 

In  this  study,  the  empirical  data  of  Snodgrass, 
et  al.  were  used,  and  were  extrapolated  to  provide  an  esti- 
mate of  the  most  likely  attenuation  at  higher  frequencies. 
The  extrapolation  of  the  attenuation  data  is  illustrated  in 
Figure  10,  a  semi-log  plot  of  the  logarithmic  attenuation 
coefficient  in  decibels  per  latitude  degree  versus  frequency 
The  empirical  data  were  extrapolated  to  the  frequency  at 
which  the  logarithmic  attenuation  function  reached  a  value 
of  1.0  dB/degree;  higher  frequencies  were  assumed  to  be 
equally  affected  by  attenuation.   Frequencies  below  0.06 
hertz  were  similarly  assumed  to  be  attenuated  by  the  same 
amount  as  the  0.06  hertz  component. 
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Figure  10:   High  Frequency  Attenuation  Function 
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Following  the  observation  of  Snodgrass,  et  al. 
that  attenuation  beyond  the  area  encompassing  the  storm  fetch 
was  negligible,  the  energy  attenuation  function  of  Snodgrass, 
et  al.  (1966)  was  applied  over  the  propagation  distance  from 
the  point-source  to  the  outer  limit  of  the  storm.   The  storm 
limit  was  taken  to  be  the  boundary  beyond  which  the  isobars 
no  longer  follow  the  circular  pattern  of  the  storm  and  thus 
can  no  longer  be  associated  with  the  counter-clockwise  flow 
around  the  low-pressure  system.   Figure  11  shows  the  location 
of  the  storm  limit  on  a  selected  sea-level  pressure  chart, 
d.   Propagation  of  Spectral  Energy 

The  energy  dissipating  effects  of  angular  spread- 
ing and  attenuation  during  propagation  of  the  swell  from  the 
point-source  on  each  successive  weather  chart  to  Monterey  is 
accounted  for  by  applying  the  spreading  and  attenuation  fac- 
tors to  the  energy-density  values  of  selected  periods  cover- 
ing the  entire  spectrum.  For  a  given  frequency  component, 
I.   ,  energy  density  propagation  loss  is  computed  by: 


[*.(«.)]'  =  [MO] 
[a.  10] 

[a  10] 


2 

-20X 

F  e 

s 


where     |_  oli'J   =  energy  density  in  deep  water  at  Monterey 

=  energy  density  at  point-source 


F$    =  angular  spreading  factor 

-  -2CLX  .       - 

e     =  attenuation  function 
By  propagating  each  of  these  dissipated  energy-density 
values  to  Monterey  at  the  appropriate  group  velocity,  the 
distribution  of  spectral  energy  with  arrival  time  is  obtained 
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for  each  point-source.   This  energy  density-time  curve  is 
referred  to  here  as  a  "propagated  spectrum".   Plots  of  the 
propagated  spectra  computed  for  deep  water  off  Monterey  for 
the  five  storms  studied  are  shown  in  Figures  12  through  16. 

In  each  figure   the  uppermost  envelope  that  can 
be  drawn  to  the  family  of  spectra  shown  is  considered  to 
represent  the  total  propagated  spectrum  for  the  actual 
swell  train  that  would  be  expected  off  Monterey.   It  may  be 
seen  in  the  case  of  Storm  1   (Figure  12) ,  that  the  spectrum 
generated  at  the  time  of  the  weather  chart  of  0000Z/16 
February  1967  had  sufficient  energy  in  all  frequencies  to 
dominate  the  other  spectra  arriving  off  Monterey.   For  the 
other  storms,  it  is  seen  that  the  swell  at  Monterey  is 
dominated  at  different  times  by  waves  arriving  from  differ- 
ent locations  along  the  storm  path. 

The  frequency  associated  with  a  particular 
energy  density  in  a  selected  propagated  spectrum  cannot  be 
read  from  Figures  12  through  16,  but  may  be  determined  from 
the  corresponding  period-time  curve  (Figures  4  through  8). 
All  spectra  shown  in  Figures  12  through  16  are  greatly 
reduced  in  energy  in  the  high  frequencies  relative  to  the 
low  frequencies  compared  to  the  spectra  present  at  the 
source-points.   The  truncated  propagated  spectra  that  may 
be  seen  in  the  figures  are  non-fully  arisen  spectra. 
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Figure  12:   Propagated  Spectra  for  Storm  1 
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Figure  13:   Propagated  Spectra  for  Storm  2 
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Figure  14:   Propagated  Spectra  for  Storm  5 
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D.   PREDICTION  OF  SWELL  CHARACTERISTICS 
1 .   Deep  Water  Period  Prediction 

Each  curve  in  Figures  4  through  8  potentially 
represents  the  period-time  variation  of  the  dominant  energy 
swell  arriving  at  Monterey  from  one  of  the  five  storms.   The 
particular  curve  or  portion  thereof  containing  the  dominant 
energy  at  a  given  moment  can  be  determined  by  referring  to 
Figures  12  through  16.   The  upper  envelope  of  the  family  of 
propagated  spectra  shown  in  each  figure  represents  the  time 
distribution  of  the  dominant  energy  waves  in  the  combined 
swell  arriving  from  all  point- sources .   The  date-time  group 
identifying  each  spectrum  refers  to  the  map  time  of  the 
point-source  from  which  the  spectrum  was  derived,  and  can 
be  used  to  relate  the  spectrum  to  a  particular  period-time 
curve  (in  Figures  4  through  8) . 

The  portions  of  those  period-time  curves  containing 
the  dominant  energy  were  identified  in  this  way  and  appear 
in  Figures  17  through  21  as  heavy  lines.   The  dots  plotted 
in  each  figure  represent  observed  period  data  and  are  dis- 
cussed in  a  later  section. 

The  reader  should  note  that  the  heavy  lines  in 
Figures  17  through  21  represent  the  predicted  distribution 
with  time  of  the  dominant  swell  periods  in  deep  water  at 
the  sensor  site.   The  effects  of  shoal-water  modifications 
to  the  arriving  swell  are  discussed  later. 
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2 .   Deep  Water  Height  Prediction 

a.  Estimation  of  Arriving  Swell  Energy 

A  moving  storm  may  be  viewed  as  the  source  of 
an  infinite  number  of  energy  spectra  generated  during  the 
storm's  lifetime.   The  family  of  energy-time  curves  for 
Monterey  shown  in  Figures  12  through  16,  then,  are  derived 
from  a  finite  number  of  selected  spectra;  specifically,  one 
spectrum  associated  with  each  weather  map  at  the  point- 
source  identified.   By  making  a  time-cut  through  one  of 
these  composites,  the  energy  densities  of  the  spectral 
components  arriving  at  Monterey  at  that  particular  time  can 
be  read  from  the  propagated  spectra.   The  period  of  the 
waves  associated  with  each  of  these  energy-densities  can 
be  read  from  the  corresponding  period-time  curve  (Figures  4 
through  8)  by  making  an  identical  time-cut.   Thus,  for  any 
selected  time,  a  plot  can  be  made  of  energy  density  versus 
frequency.   A  continuous  spectrum  drawn  through  the  points 
of  such  a  plot  represents  the  predicted  energy  spectrum  of 
the  swell  in  deep. water  at  Monterey  at  that  time.   Figure  22 
illustrates  these  predicted  spectra  for  a  series  of  arrival 
times  for  Swell  Train  3.   The  area  under  each  spectrum  is 
proportional  to  the  predicted  energy  in  the  arriving  swell. 

b.  Calculation  of  Swell  Heights  from  Energy  Estimates 
The  total  energy  in  a  sea  is  found  by  integrat- 
ing the  energy  spectrum  over  the  full  range  of  frequencies 
from  0  to  oo  .   This  energy  value,  if  the  approximate  spec- 
tral form  is  known,  can  also  be  related  to  the  statistical 
distribution  of  wave  heights  in  the  sea.   Longuet-Higgins 
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(1952)  has  shown  that  for  narrow  band  spectra  (i.e.,  swell) 
the  significant  wave  height,  the  average  of  the  highest 
third  of  the  waves  in  a  wave  record,  is  given  by: 

(4) 


H,  =  2.83  V  E 


where     H  =  significant  wave  height 

E  =  total  energy  in  the  spectrum 
If  the  predicted  swell  spectra  are  sufficiently  narrow 
banded  in  the  sense  implied  by  Longuet-Higgins  (1952)  , 
this  statistical  relationship  provides  a  means  for  computing 
the  swell  height  distribution  with  time  from  the  predicted 
energy-time  distribution. 

3 .   Refraction  and  Shoaling  Modifications 
a.   Shoal  Water  Height  Prediction 

The  wave  sensor  used  to  obtain  the  observed 
swell  height  at  Monterey  for  the  five  storms  studied  was 
located  in  shoal  water;  accordingly,  the  effects  of  refrac- 
tion and  shoaling  had  to  be  considered  in  order  to  compare 
observed  heights  with  heights  derived  from  the  computed 
swell  energy. 

Since  refraction  effects  are  dependent  on  the 
direction  of  swell  arrival,  and  directional  information  was 
not  recorded,  two  procedural  alternatives  presented  them- 
selves:  the  arrival  direction  of  the  observed  swell  could 
be  estimated,  and  the  observed  heights  adjusted  to  remove 
the  effects  of  refraction  and  shoaling;  or  the  predicted 
heights  could  be  modified  for  these  effects  and  compared 
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with  the  shoal  water  data  as  observed  at  the  sensor  site. 
The  latter  course  of  action  was  chosen  because  an  arrival 
azimuth  for  the  swell  train  propagated  from  each  point- 
source  is  a  function  of  the  great  circle  path  between  the 
point-source  and  Monterey. 

To  illustrate  the  procedure  used  to  calculate 
the  height  distribution  in  the  swell  at  the  wave-gage  site, 
Storm  3  will  be  used  as  an  example  (Figure  22) .   The 
individual  energy-density  values  determined  from  the  deep 
water  spectral  time-cuts,  indicated  by  the  dots  in  the 
figure,  are  modified  for  the  effects  of  refraction  and 
shoaling  using  the  relation: 

[asi.,)]2  =  [a.(«,i]Vk,  (    (5) 

where      I".  /  .  il  =  shoal-water  energy  density  associated 

with  frequency  f. 

Ao(f.)  =  deep-water  energy  density 

Kr  =  refraction  coefficient 

Ks  =  shoaling  coefficient  (water  depth  30 

feet) 

Refraction  coefficients  for  the  sensor  site  are  shown  in 

Figure  23;  the  data  were  compiled  for  the  sensor  site  from 

a  number  of  manually  drawn  refraction  diagrams  prepared  as 

a  laboratory  exercise  by  former  students  at  the  Naval 

Postgraduate  School. 

The  energy  in  the  swell  arriving  at  the  sensor 

site  at  a  selected  time  is  then  approximated  using  the  modi 

fied  energy  density  values  in  the  step-wise  integration 

scheme : 
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n  — i 


n  — i 

T    (6) 


k(ol 


1=1 

where   E  (f)   =  total  shoal-water  energy  in  swell  arriving  at 

time  t 
=  modified  energy-density  value  associated  with 
frequency  f. 
n   =  number  of  energy-density  values  used  in  the 

approximation 
Al   =  frequency  band  between  I   and  f 
Figure  24  illustrates  the  shoal  water  energy  spectra  for 
Swell  Train  3  at  selected  times. 

Predicted  values  of  significant  swell  height  in 
shoal  water  may  then  be  computed  from  these  shoal-water 
energy  computations  using  Equation  4.   Predictions  of  signi- 
ficant wave  height  for  the  swell  arriving  at  the  wave  gage 
site  from  Storm  3  are  shown  in  Figure  26. 
b.   Shoal  Water  Period  Prediction 

As  a  result  of  refraction  and  shoaling,  the 
energy  density  associated  with  a  given  period  may  be  expected 
to  differ  at  a  shoal  water  site  from  the  deep  water  value. 
Because  these  processes  are  frequency  dependent,  a  shift  of 
the  dominant  energy  may  occur  from  one  period- time  curve  to 
another  as  the  swell  passes  from  deep  to  shoal  water. 

A  prediction  of  the  dominant  periods  in  shoal 
water  is  obtained  by  plotting  the  periods  of  maximum  energy 
density  in  the  shoal  water  spectra  as  a  function  of  arrival 
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time.   Figure  28  illustrates  the  predicted  shoal  water 
period- time  curve  for  Swell  Train  3,  as  determined  from  the 
shoal  water  energy  spectra  shown  in  Figure  24. 
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III.   COMPARISON  OF  PREDICTED  AND  OBSERVED  SWELL 

A.   OBSERVED  WAVE  DATA 

The  observed  wave  data  used  in  this  study  for  comparison 
with  the  computed  wave  data  were  obtained  from  analysis  of 
strip-chart  recordings  from  two  bottom-mounted  pressure- type 
wave  sensors  located  in  approximately  30  feet  of  water.   The 
strip-chart  records  from  both  sensors  were  manually  analyzed 
by  Professor  Warren  C.  Thompson.   Wave  periods  were  obtained 
using  the  wave-group  method,  which  for  distant  swell  has 
been  shown  to  yield  periods  equivalent  to  the  periods  of 
maximum  energy  density  obtained  by  spectral  analysis 
(Thompson,  1973) .   The  wave  records  exhibiting  the  swell 
originating  from  Storms  1  through  4  were  obtained  from  the 
Monterey  sensor  located  at  Del  Monte  Beach.   The  swell 
periods  from  Storm  5  were  obtained  from  a  sensor  at  Stinson 
Beach,  California.   The  Stinson  Beach  sensor  is  located 
about  75  nautical,  miles  closer  to  the  swell  source  than  the 
Monterey  sensor  and  approximately  on  the  same  great  circle 
trajectory.   Since  the  predicted  periods  were  computed  for 
Monterey,  the  arrival  times  of  the  observed  periods 
recorded  at  Stinson  Beach  were  adjusted  to  correspond  to 
arrival  times  at  the  Del  Monte  Beach  site. 

Significant  wave  heights  at  the  Monterey  gage  were 
determined  by  adjusting  the  apparent  heights  recorded  by 
the  pressure  gage  for  the  effects  of  hydrodynamic  damping 
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due  to  water  depth  using  Wiegel's  Tables  (Wicgcl,  1954). 
The  recorded  swell  was  presumed  to  be  narrow  banded  so  that 

a  single  damping  factor  could  be  applied  with  good  approxi- 
mation.  The  pressure-response  factor  of  the  wave  gage  was 
greater  than  0.988  over  the  period  range  from  5  to  20 
seconds.   The  significant  wave  heights  computed  from  the 
pressure  gage  data  are  referred  to  as  observed  heights. 

B.   SWELL  HEIGHT 

Significant  swell  heights  at  the  Monterey  gage  site 
were  predicted  for  the  swell  arriving  from  two  of  the  storms 
studied  (Storm  1  and  Storm  3) .   Comparisons  of  predicted  and 
observed  signficant  heights  are  shown  in  Figures  25  and  26. 
It  may  be  seen  that  in  both  sets  of  swell  the  observed 
heights  are  greater  than  the  predicted  heights;  in  Swell 
Train  1  the  observed  heights  exceed  those  predicted  by  a 
factor  of  about  2.6,  whereas  in  Swell  Train  3  the  factor  is 
about  5.   In  the  case  of  Swell  Train  1,  the  observed  peak 
height  occurred  approximately  ten  hours  earlier  than  pre- 
dicted.  The  peak  height  in  Swell  Train  3  occurred  approxi- 
mately at  the  time  predicted. 

The  difference  in  magnitude  between  the  predicted  and 
observed  heights  most  likely  can  be  attributed  to  one  or 
more  of  the  procedures  used  to  apply  angular  spreading, 
attenuation,  and  shoaling  and  refraction.   The  choice  of  a 
12-millibar  interval  for  defining  the  fetch  width  is  largely 
arbitrary  and  may  be  too  conservative.   The  way  in  which 
the  storm  limit  is  defined  may  have  introduced  excessive 
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attenuation  estimates.   Poor  refraction  data  also  may  have 
introduced  a  significant  error  in  the  calculations  of  the 
shoal-water  swell  energy. 

C.   DOMINANT  SWELL  PERIOD 

The  distribution  with  time  of  the  dominant  swell  period 
at  the  wave  gage  site  also  was  computed  for  Swell  Trains  1 
and  3.   The  results  are  shown  in  Figures  27  and  28  by  open 
circles;  the  observed  periods  are  represented  by  dots.   The 
predicted  and  observed  periods  at  the  sensor  site  are  super 
imposed  on  the  family  of  potential  period-time  curves  com- 
puted for  Monterey  as  shown  in  Figures  4  and  6.   As  stated 
earlier,  the  heavy  portions  of  the  period-time  curves 
represent  the  predicted  dominant  periods  in  deep  water  at 
Monterey . 

In  Figure  27  it  may  be  observed,  from  the  predicted 
period  distribution  with  time  in  deep  water  and  at  the  wave 
gage  site,  that  the  shoaling  and  refraction  modifications 
to  the  propagated  spectra  resulted  in  a  shift  of  the  domin- 
ant energy  from  one  period-time  curve  to  another  in  the 
period  band  from  18  to  21  seconds.   Figure  28  exhibits  a 
similar  shift;  however,  in  this  case  the  predicted  shift  is 
more  consistent  with  the  trend  of  the  observed  periods.   It 
may  be  concluded  that  the  dominant  period  can  shift  due  to 
shoaling  and  refraction,  but  that  the  shift  is  not  large. 
Accordingly,  the  predicted  period-time  distributions  in 
deep  water  are  used  for  comparison  with  the  observed  peri- 
ods at  the  wave  gage  for  Storms  2,  4,  and  5. 
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Comparison  of  the  period- time  distribution  of  the  pre- 
dicted dominant  swell  (heavy  lines)  with  the  mean  trend  of 
the  observed  periods  (dots)  for  each  of  the  storms  (Figures 
17  through  21)  shows  close  agreement.   For  the  five  storms 
studied,  the  mean  difference  between  the  predicted  and 
observed  periods  does  not  exceed  about  one  second.   It  may 
be  concluded  that  the  method  produces  accurate  period 
forecasts . 

Differences  in  the  swell-generating  characteristics 
among  the  five  storms  may  be  noted  from  examination  of  the 
predicted  periods  (heavy  lines)  and  the  observed  periods 
(dots)  shown  in  Figures  17  through  21.   In  Storm  1  (Figure 
17),  the  period-time  curve  propagated  from  a  single  point- 
source  (on  the  weather  chart  of  0000Z/16  February  1967) 
contains  the  dominant  energy  over  the  entire  duration  of  the 
observed  swell  at  Monterey;  this  is  verified  by  the  close 
correspondence  of  the  predicted  deep-water  period  distribu- 
tion with  the  observed  periods.   In  Storm  4  (Figure  20), 
the  distribution  of  observed  periods  also  indicates  effec- 
tive generation  at  a  single  point-source;  however,  predic- 
tion called  for  the  longest  period  components  of  the  swell 
to  be  generated  at  the  point-source  identified  on  the 
succeeding  weather  map,  and  to  arrive  sooner.   It  is  evident 
from  an  examination  of  the  plot  of  propagated  spectra 
derived  from  this  storm  (Figure  15) ,  that  determination  of 
which  period- time  curve  is  dominant  at  a  given  time  is  in 
this  case  particularly  dependent  upon  the  accuracy  of  the 
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modifications  to  the  sea  spectrum  at  each  point-source 

due  to  spectral  cutoff,  angular  spreading,  and  attenuation. 

By  contrast,  in  Storm  3  (Figure  19)  there  is  a  notice- 
able shift  of  the  predicted  dominant  periods  as  well  as  the 
observed  periods  from  one  period-time  curve  to  another  such 
that  the  dominant  energy  arriving  at  different  times  is 
propagated  from  several  different  point-sources.   In  Storm  2 
(Figure  18) ,  a  similar  shift  in  the  dominant  periods  across 
the  dispersion  curves  may  be  noted.   In  this  storm,  the 
period-time  curves  are  nearly  coincident,  which  is  due  to 
the  fact  that  the  speed  of  the  point-source  toward  Monterey 
and  the  group  velocity  of  the  peak  energy  component  in  each 
source  spectrum  were  closely  coincident.   The  wide  scatter 
of  periods  recorded  after  0000Z/25  February  is  attributed 
to  a  wave  train  from  another  disturbance. 

D.   REFINED  PERIOD  PREDICTION 

In  order  to  further  investigate  the  apparent  shift  of 
the  dominant  energy  from  one  period- time  curve  to  another, 
which  is  particularly  apparent  in  the  case  of  Storm  3 
(Figure  19) ,  the  synoptic  parameters  derived  from  each  sur- 
face pressure  chart  (i.e.,  surface  wind  speed,  storm  speed, 
and  distance  from  the  point-source  to  Monterey)  were 
linearly  interpolated  at  two-hour  increments.   The  period- 
time  distribution  of  the  dominant  swell  at  Monterey  obtained 
from  this  refined  data  base  is  illustrated  in  Figure  29  by 
the  heavy  lines.   The  associated  propagated  energy  spectra 
are  illustrated  in  Figure  30.   The  reader  may  note  that  the 
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time  sequence  of  the  dominant  segments  of  the  period- time 
curves  reveals  the  generation  history  of  the  arriving  swell 
train.   In  the  period  band  from  18  to  23  seconds,  the  domin- 
ant energy  was  generated  in  the  later  part  of  the  storm 
history  (1000Z/13  November  through  2000Z/13  November);  the 
period  band  from  13  to  18  seconds  was  generated  early  in  the 
storm  history  (1200Z/12  November  through  1800Z/12  November) . 

The  general  trend  of  the  dominant  periods  in  the  pre- 
dicted swell  may  be  seen  to  approximate  the  trend  of  the 
observed  periods;  however,  the  correlation  does  not  appear 
to  have  been  significantly  improved  by  refinement  of  the 
six-hourly  synoptic  data. 
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Figure  30:   Refined  Propagated  Spectra  from  Storm  3 
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IV.   CONCLUSIONS 

In  spite  of  the  poor  height  predictions  obtained  for 
the  two  swell  trains  at  the  Monterey  wave  gage  site,  the 
multiple  point-source  method  is  sufficiently  promising  with 
regard  to  the  time  of  arrival  of  the  maximum  swell  height 
and  the  accuracy  of  the  dominant  periods  obtained  that 
additional  research  effort  should  be  expended  to  test  the 
method  further,  and  to  tune  the  procedures  involved  so  as 
to  give  better  height  predictions.   The  latter  may  be 
accomplished  by  modifying  some  of  the  procedures  within 
realistic  limits,  particularly  the  procedures  for  angular 
spreading  and  attenuation.   The  predicted  heights  should 
be  compared  with  observed  wave  data  of  better  quality  than 
were  available  in  this  study,  preferably  recorded  in  deep 
water  so  that  the  complications  of  shoaJLing  and  refraction 
can  be  eliminated. 
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